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Asymptotically flat, time-symmetric, axially symmetric 
and conformally flat initial data for vacuum general relativity 
are studied numerically on R 3 with the interior of a standard 
torus cut out. By the choice of boundary condition the torus is 
marginally outer trapped, and thus a surface of minimal area. 
Apart from pure scaling the standard tori are parameterized 
by a radius a £ [0, 1], where a = corresponds to the limit 
where the boundary torus degenerates to a circle and a = 1 to 
a torus that touches the axis of symmetry. Noting that these 
tori are the orbits of a U(l) x U(l) conformal isometry allows 
for a simple scheme to solve the constraint, involving numer- 
ical solution of only ordinary differential equations. Regular 
solutions can be constructed for a €]ao ~ 0.14405, 1[. The 
tori are unstable minimal surfaces (i.e. only saddle points of 
the area functional) and thus can not be apparent horizons, 
but they are always surrounded by an apparent horizon of 
spherical topology, which is analyzed in the context of the 
hoop conjecture and isoperimetric inequality for black holes. 

PACS numbers: 4. 20. Ex, 2.40.Ky, 4. 25. Dm 



I. INTRODUCTION 

The aim of this paper is to construct a class of asymp- 
totically flat solutions to the vacuum constraint equations 
of general relativity that are somewhat unusual, namely 
they contain a marginally outer trapped torus. In the 
case of time symmetry discussed here, a marginally outer 
trapped surface has vanishing mean curvature and is a 
surface of minimal area (see Ref. Q and references cited 
therein). The minimal tori are all unstable, that is, they 
are only saddle points of the area functional. In a re- 
cent discussion of minimal area tori and their connection 
with topological censorship Galloway fl proved, that an 
unstable minimal torus locally can be deformed into an 
outer trapped surface. The outer boundary of the re- 
gion containing trapped surfaces - the apparent horizon 
- is a stable marginally outer trapped surface of spherical 
topology (see Sec. |v|). 

Motivation for this work comes from two sources: while 
the question of the existence of a regular solution to the 
constraints is of physical and mathematical significance 
in itself, analysis of given initial data with respect to the 
existence and properties of minimal surfaces and appar- 
ent horizons has received some interest in the discussion 



on cosmic censorship. 

This work deals with the simplest situation where a 
toroidal minimal surface can be present, namely that of 
time-symmetry, axial symmetry and conformal flatness. 
The assumptions of time-symmetry and conformal flat- 
ness allow to use a conformal symmetry to turn the con- 
straint from a partial differential equation (PDE) to a 
set of uncoupled ordinary differential equations (ODEs). 
Generalizations of the scheme that relax the restriction to 
conformal flat nes s and axial symmetry will be discussed 
below in Sec. 



VI 



The physical data are given by an asymptotically fiat 
metric g defined on a manifold M that is taken to be R 3 
with a region inside a torus dM cut out. The physical 
metric g is required to satisfy the vacuum Hamiltonian 
constraint of vanishing scalar curvature R g on M plus 
the boundary condition of vanishing mean curvature p/2 
(p is the trace of the extrinsic curvature) at dM, making 
it a minimal surface: 



Rg = on M, 



on dM. 



p = Q 



(1) 
(2) 



The standard way to solve these equations is the con- 
formal approach ||. Restricting to conformal flatness 
we select the flat metric S = dr + dz + r 2 dip on M 
as our base metric and construct the physical metric as 
g — ip 4 S, which reduces Eq. (|l|) to an elliptic PDE for 
the conformal factor ip. 

For notational convenience we define the linear differ- 
ential operators 



A s +\r 



li- 



the conformal Laplacian of a metric <?, and 

d 1 



B„ 



dn 



(3) 



(4) 



where d/dn is the normal derivative and p is the trace 
of the extrinsic curvature (or twice the mean curvature) 
induced by the 3-metric jona hypersurface. Using the 
definitions of L g and B g the constraint equation ([!]) with 
boundary condition (^|) turns into the problem of finding 
a conformal factor ip, such that 
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Lgt/j = - A s ip = on M, 
B s ip = on dM, 
ip > 0, V 



1 as r + z 



(5) 

(6) 
ex.. (7) 



We are left with the freedom to fix the boundary sur- 
face, taking it to be a standard torus in flat space, defined 
in cylindrical coordinates as a solution of 



z 2 + (r - Af 



a 2 . 



(8) 



Since a finite overall scaling of the coordinates bears no 
physical significance, we can fix A = 1 for the boundary 
tori considered, which leaves us with the one parameter 
family < a < 1. The extremal case a = corresponds 
to a circle of unit radius in the z = plane, the other 
extreme, a = 1, defines a torus which touches its axis of 
symmetry. 

Using the method described below it is possible to con- 
tinue the data smoothly inside the minimal tori to some 
finite extent, assuming vacuum and conformal flatness in- 
side, but in general this assumption will lead to a nodal 
surface, and thus to a curvature singularity, if one pro- 
ceeds too far. The obvious question whether it is possible 
to get a regular interior and nonnegative energy density 
of the matter is not dealt with here. 

Regular solutions to the constraint equation (||) with 
boundary conditions (|L [?]) exist for a e]a ~ 0.14405, 1[ 
(all numerical results are rounded to the last digit), where 
the parameter a labels the boundary tori as defined above 
in Eq. (||). In the limit a — > ao the mass blows up, and 
the Sobolev quotient (an extension of the definition of the 
Yamabe number to manifolds with boundary) vanishes. 
Also the apparent horizon moves outward to infinity and 
becomes spherical in the limit. This behavior is famil- 
iar from the case when the (usual) Yamabe number ap- 
proaches zero ||. If a is further decreased, the Sobolev 
quotient becomes negative, the mass diverges to minus 
infinity as a / a , and the conformal factor develops a 
nodal surface. In the other extremal case, when a — ► 1, 
the mass stays finite, but a curvature singularity develops 
at the center, where the toroidal minimal surface touches 
the axis of symmetry. In this last case the solution can 
be given explicitly. 

The toroidal minimal surfaces are all unstable, and 
part of their spectrum of perturbations was computed. 

The ADM mass can be cast into an explicitly positive 
form in a way similar to Brill's result Q] , additionally it 
can be bounded from below. 

Note that the initial data discussed here are perfectly 
viable for evolution, since an apparent horizon boundary 
condition f| can be applied at the topologically spheri- 
cal apparent horizon that surrounds the boundary torus. 
The minimal torus could then be viewed as a model of a 
bulk of matter or gravitational radiation that has enough 
energy to confine all information about it inside of the 
horizon. This comment also holds in the case a — * 1, 
when a curvature singularity appears at the center. 



By using conformal symmetries, the PDE can be de- 
composed into an infinite set of uncoupled linear ODEs 
that were solved numerically. A convenient cutoff of the 
infinite series yields a highly accurate solution for the 
constraint. The existence of solutions can be decided by 
solving a single linear ODE, together with the high ac- 
curacy obtained for the actual solution, this enables a 
secure grip on the two interesting limiting cases. 

The key technical input is the utilization of a U(l) x 
U(l) conformal isometry that leaves the boundary torus 
invariant. In order to benefit from the conformal isome- 
try the conformal method is reformulated for a compact 
base manifold M, taken to be the one-point compactifi- 
cation of M, that is S 3 with the interior of a torus cut 
out. The conformal ansatz g = G 4 g, with g the standard 
metric on S 3 , results in an elliptic PDE similar to (H), 



L g G 



Aire S a 



on M, 



(9) 



but with a distributional source term: (5a is the Dirac 
delta distribution concentrated at the point A on M, 
which corresponds to asymptotic infinity of M, and c 
is a positive constant that can be chosen arbitrarily for 
overall scaling. G is thus a Green function for c~ x L g with 
boundary condition 



B g G = 



on dM. 



(10) 



The existence of a positive Green function G is deter- 
mined by the sign of Ai , the first eigenvalue of the con- 
formal Laplacian with the present boundary condition, 
or equivalently by the sign of the Sobolev quotient Q(g), 
a conformal invariant defined as 



Q(g) ■-- 



inf 



Jm(IWI 2 + \RgV 2 )dv + \ J dM PV 2 dv 



(11) 



which is a generalization of the Yamabe functional || 
(defined for manifolds without boundary) to manifolds 
with boundary |t]] and has the same sign as Ai. 

It can be shown in a manner equivalent to the no- 
boundary case dealt with in || , that a positive G exists 
iff Ai is positive. Furthermore, when G > 0, the positive 
mass theorem holds for g = G 4 g 0. 

The rest of this paper is organized as follows: In the 
next section (||) conformal compactification via inverse 
stereographic projection from R 3 to S 3 and toroidal coor- 
dinates thereon will be discussed. In Sec. Ill the method 



to solve the constraint numerically using toroidal coor- 
dinates will be developed, including numerical solution 
for the first eigenvalue of the conformal Laplacian. Sec- 
tion IV deals with the instability of the minimal tori by 
considering normal variations of the surfaces. Apparent 
horizons are located and analyzed with respect to the 
isoperimetric inequality for black holes [ fl0[ and the hoop 
conjecture in Sec. [v| A discussion of the results, 
with emphasis on the limiting cases, is given in section 
VI. 
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II. COMPACTIFICATION AND CONFORMAL 
SYMMETRY 

For every choice of boundary torus £ 7 , that is a torus 
described by Eq. (||) with A = 1 and parameterized by 
7 := \J\ — a 2 < 1, we define a conformal rescaling from 
the flat metric 5 on R 3 to a metric (7(7) defined on S* 3 by 



r, (12) 



3(7) = £l(r,z; 7) 2 (5, »(r,z; 7 ) = 

7 Z + 7' z -r 

which corresponds to an inverse stereographic projection. 
The notation f(r, z\ 7) is chosen to denote dependence of 
a function / on both the spatial coordinates and the pa- 
rameter 7, frequently this will be shortened to writing 
7(7). With the chosen normalization 5(7) is the stan- 
dard metric on the 3-sphcrc of radius 1/7, which is con- 
veniently written as 



fl(7) 



dr 



+ (l-ff^df+fV 



1 — <y j" 



(13) 



in toroidal coordinates (f , z, ip) taking ranges < f < 
1/7: < "fz, ip < 2tt. A discussion of the construction 
of this coordinate system on S* 3 of unit radius is given 
in |[2|. The coordinate f as a function on the (r, z)-half 
plane can be read off from Eqs. ( |l2] ) and (|l3|) yielding 

r(r, z; 7) = f2(r, z; 7) r. 

The relation between the coordinate z and the original 
cylindrical coordinates is 



cos (jz) 



1 -7 2 »(r,z; 7 ) 
V 7 ! -i 2 f(r,z\'fj 1 



sign(z) = sign(z) 



The point A on S 3 , corresponding to asymptotic infinity 
of R 3 , is given by f = z = 0. 

The vector fields (d/dp) a , (d/dz) a form a pair of com- 
muting, orthogonal and hypersurface-orthogonal Killing 
vector fields of 17(7) spanning the surfaces of constant f. 
For f ^ 0,1/7 these are fiat tori of constant mean cur- 
vature changing sign at the so called Clifford torus for 
which f = I/7V2. Going back to the flat metric and 
cylindrical coordinates, a torus f = const, is mapped to 
a standard torus as defined by Eq. (||) with parame- 
ters A = 1/f, a = yl — 7 2 f 2 /f, in particular the torus 
f = 1 is mapped to the boundary torus £ 7 . The sets 
f = (respectively f = I/7) are linked great circles on 
5 3 , corresponding to the z-axis (respectively the circle 
z = 0, r = 7) after stereographic projection. 

The totally geodesic surfaces z = const, (ip — const.) 
are metric hemispheres intersecting at f = I/7 (f = 0) 
with {jz — a}L){jz — a+n} ({tp — (3}U{p = f3+ir}) be- 
ing smoothly embedded 5 2 's. Under stereographic pro- 
jection, the surfaces 72 = const, do not change topology 
except for {72: = 0}U{7Z = ir} which gets decompactified 
into {z = 0, r > 7} U {z = 0, r < 7}, i.e. the equatorial 



plane z — 0. In flat space {72: = a} U {72 = a + tt} 
is a sphere centered at r = 0, z = 7Cot7z with radius 
yV (cot 2 jz- 1) + 2. 

Note that for 7 = the radius of S 3 becomes infinite, 
corresponding to i? 3 , and the coordinates f, z become 
inverted cylindrical coordinates thereon: 



2r 



2z 



r 2 + z A 



The tori f = const., corresponding to A = a = 1/f, all 
touch the z-axis in this limiting case. 

We have thus constructed a family of foliations of R 3 
by standard tori, with each element defined by selection 
of a boundary torus S 7 that is contained in the foliation. 
For every foliation labeled by 7 we defined in Eq. ( |l2| ) a 
conformal rescaling with compactifying factor f2(r, z;j). 
Since f2 is independent of (p, (d/dtp) a is a Killing vector 
of g, but (d/dz) a is only a conformal Killing vector of g. 

The reason for the chosen scaling of S 3 as opposed to 
working on S 3 of fixed radius can bee seen when rescaling 
to the standard metric g on S 3 of unit radius (with the 
notation of |l3|): 



9 ■= i 2 g(i) 



dp 2 



+ (i-p 2 )d x 2 + P 2 dip 2 



1-p 2 

using rescaled toroidal coordinates (p, x, p) 

p(r,z;j) :=jf(r,z;j), x( r , z ;i) -=iz. 

taking ranges < p < 1, < x,<P < 27r. The boundary 
torus S 7 that was defined by f = 1 before is now specified 
by P — li so that the region of physical interest, < p < 
7 shrinks to zero volume in the limit 7 — > 0. The problem 
of finding a conformal factor solving the rescaled equation 
LgG = 47tv / 2<5a then becomes ill-defined, the conformal 
factor G (and thus the mass of G 4 g) blows up for 7^0. 

In order to formulate the constraint equation explicitly 
we need to write down the operators L g and B g in our 
toroidal coordinate system. The extrinsic curvature of a 
torus of constant f with respect to the metric g a b is (the 
semicolon denotes the covariant derivative as usual) 



p ab = r\J\- 7 2 r 2 (p-aP-b ~ l 2 z. a z. b ) 
At dM (f = 1) this yields 

!- 2 7 2 

P = Pa ~ 



(14) 



(15) 



which vanishes at the Clifford torus given by j c = 1/ y/2, 
where dM is maximal. The extrinsic curvature p ao with 
respect to the physical metric g is 



Pab = G 2 {pab + 2g ab n c G- c ) , 



(16) 



where n c is the unit normal vector field on dM. Note 
that since p ab has Lorentz signature, in particular it is 
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not proportional to g ([T||), by the conformal rescaling 
formula for the extrinsic curvature ( |l6|) it is not possible 
to get p ao — by a conformal transformation. As a 
consequence the minimal tori cannot be stable, since in 
that case they would have to be totally geodesic (p a b = 0) 
- (see |). 

Finally we can write down the operators L 7 := L 9i 
and -B-v := B g - simplifying their labels - defined in 



(17) 



Eqs. (g, [tC]) explicitly: 



B 7 = VI 




-Zv-y 



(1 



2-2, d 2 
d 2 



1 



1 - 37 2 f 2 9 
1 d 2 3 , 



1 — 7 2 r 2 <9z 2 



<9(,2 2 



(18) 



V- = i + 



2Vr 2 + z 2 



o(si). 



The limiting cases 7 = 0, 1 can be considered as singu- 
lar situations, and the question arises how this will effect 
the physical metric 5(7), in particular its existence. For 
7 < 7 C we have p positive (see Eq. ( n5 l )) w hence is the 
Sobolev quotient Q(g) defined in Eq. Qll[) - a positive 
Green function thus exists and the mass of the physical 
metric 3(7) is positive. On the other hand p is negative 
for 7 > 7c, so that a change of sign becomes possible for 
Q{g) at some 7, 7 C < 7 < 1. This indeed happens at the 
numerically calculated value 70 = 0.98957, correspond- 
ing to a = 0.14405. As one expects, this case is analogous 
to the Yamabe number going to zero in the usually con- 
sidered case of manifolds without boundary. The mass 
diverges to plus infinity as 7 approaches 70 from below, 
and to minus infinity when taking the limit from above. 

For 7^0 the metric g becomes the flat metric on the 
infinite cylinder < r < 1, —00 < z < 00, < <p < 2ir: 



III. SOLUTION OF THE CONSTRAINT 
EQUATION 

A. Formulation 

From the choices made above our physical data are 
determined by the single parameter 7, which labels the 
position of the boundary torus. This defines the phys- 
ical metric (7(7) = G(7) 4 g(7), where G("f) satisfies the 
equations 



L 7 G(7) =4irV25 A on M, 
£ 7 G( 7 ) = on dM 



(19) 
(20) 



with operators L 7 and B 1 given by the expressions 
(0>@)) the boundary torus is located at the fixed co- 
ordinate value r = 1 and the scale factor is chosen as 
c = \/2. The rescaling to flat space defined in Eq. (JjT 
corresponds to a Green function Go, 



G (7) = O( 7 )- 1 / 2 =7 1-v/T 



7 2 r 2 cos "fz 



-1/2 



where Go (7) satisfies Eq. ( |H| ) but with the boundary 
condition ( f20[ ) replaced by regularity on 5 3 \{A}. The 
asymptotic behavior of G near A is known H to be 



G = G + i ^= + O(G - 1 ), 



(21) 



where m is the ADM-mass of the physical metric g. By 
the positive mass theorem [|l3| m is positive provided Ai 
is positive. 

The relation to the original formulation of the con- 
straint problem in Eqs. ((^-Q) is as follows: 



G 4 g = i> 4 5 = t/j 4 n- 2 g ip = fT 1/2 G, 



g{0) = df 2 +dz 2 + f 2 d<p 2 . 

The point A, representing physical infinity, is at f = z = 
0, the unphysical infinity, z — ±00 corresponds to the 
physical center r = z = 0. The constraint Eq. (JsJ) re- 
duces to the cylindrically symmetric Poisson equation for 
a point source 



L G(0) 



d 2 Id d 2 
df 2 f df dz 2 



G(0) = 4ttV2Sa, 



with boundary condition 



(22) 



(23) 



The explicit limiting solution is derived in Appendix |A| 



as 



V V 2 + z 2 



cos(wz) — - — — — — — l Q {LLir)duj. 

As is discussed in Appendix [A], G(f, z; 0) is a regular 
positive function for < f < 1, —00 < z < 00, f 2 + z 2 > 
0. The ADM-mass of g(0) = G(0) 4 5 (0) is m = 3.877, in 
particular it is finite. For \z\ large G decays exponentially, 

G(f, z; 0) oc exp(— a\z |), a = 0.6856, 

it is thus zero at the physical origin r = z — 0. There the 
Kretschmann invariant diverges as R a bcdR abcd [g(0)) oc 
exp(8a|z|), the physical metric thus has a curvature sin- 
gularity at the origin, where the minimal torus touches 
the axis of symmetry. 
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B. Existence of Solutions 



As stated above, the existence of a solution to the con- 
straint, that is a positive Green function for the elliptic 
operator L 7 with boundary condition given by Eq. (10), 
is equivalent to positivity of the lowest eigenvalue Ai of 
L 7 . Calling the eigenfunctions 6 we obtain the eigenvalue 
problem 

Lj9 — X9, B~f0\f = i = 0, 

which is separable and the solutions can be written as 

6\ mn (f, <p, z) = $\ mn {r) sa\(nip + a) sin(m7z + (3), 

where a and (3 are arbitrary angles and we are left to 
solve an equation for i3, 



1 



where L f7 is the 'radial' part of the conformal Laplacian 
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FIG. 1. The lowest eigenvalue Ai of the conformal Lapla- 
cian is plotted as a function of 7, the insert shows that Ai 
gets negative for 7 > 70. The value for 7 = is marked by 
an x. 



^ = -(l-7 2 f 2 )|s 



l- 3 7 2 f 2 d 3 2 



Since the first eigenfunction is non-degenerate jljl ] and 
has therefore to depend trivially on the angular coordi- 
nates (p and z, we may set m = n = when looking for 
the lowest eigenvalue Ai : 



Ai? 



A00- 



(25) 



The existence question has thus been reduced to an 
ODE of Sturm-Liouville type which can be solved by 
standard methods. In this case routine d02kef from the 
NAG Fortran Library Jll| was chosen to calculate eigen- 
values and eigenfunctions. 

The lowest eigenvalue is found to be positive for 7 < 
70 = 0.98957, and negative for 7 g reater, see figure (|l|) 
and table (I). For 7 = Eq. ( |25| ) reduces to Bessel's 
equation: 



fV + ftf' + Af 2 t? = 0, 



tf' + ^\r=l = 0. 



The regular solution is (apart from arbitrary scaling) -d = 
Jo(VXr) (J is the Bessel function of order 0) where A is 
chosen to satisfy the boundary condition 

Jo(V\r) + ~J (V\f)\f=i =0, 
which yields A x = 0.4700. 



C. Transforming the Constraint to ODEs 

Since the singular behavior of G(f,z;"f) at A makes a 
direct numerical solution for the Green function incon- 
venient, the first step in solving the system (^, |l0|) is to 
regularize equation (0) . We split G into the singular part 
Go and a regular part (f> by 

G(7) = G (7) + #7). (26) 
The asymptotic behavior of 4> near A is 

from the asymptotic behavior of the Green function (|2l]) , 
m is the ADM-mass of the physical metric g. Calculation 
of the mass is thus as simple as evaluating the numerical 
solution (p at the point A. Insertion of ( po] ) into (g, |u]) 
yields an elliptic problem for <fi: 




l-2 7 2 
2(1 -7 2 ) 



6(2,7). 



Note that the boundary condition becomes singular for 
7 ->l. 

The above partial differential equation in two dimen- 
sions can be turned into a system of a countable set of 
uncoupled ODEs by decomposing <f> and b into a cosine 
series (the sine terms of the Fourier decomposition are 
zero due to the symmetry with respect to reflection at 
the equator): 
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= ^2 <j> n (r) cos(n~fz), b = ^ b ^( r ) cos(n-yz). 

(27) 



n=0 



This gives a boundary value ODE problem for every n 
(prime denoting differentiation with respect to f): 



L r 



2 — <y2^,2 

27 2 







4(1-7*) 



r=i 



,o=0. 



(28) 

(29) 
(30) 



The last condition ensures regularity on f = 0, the 
z— axis. By integrating the ODEs beyond f = 1 it is pos- 
sible to continue the solution smoothly inside the minimal 
tori to some finite extent, but this will eventually lead to 
a nodal surface, and thus to a curvature singularity. 

Note that, for n = 0, Eq. ( p8| ) corresponds to Eq. (25), 
and correspondingly only <j)Q blows up for 70, the higher 
order terms pass through Ai = smoothly. 

The equations are solved by a shooting and match- 
ing method (NAG routine d02agf for n = 1, . . . , N, 
where N is suitably chosen to make the truncation error 
for ip reasonably small. 



with cutoff-parameter N one may ask how well the 
Hamiltonian constraint is satisfied. Due to numerical 
inaccuracies and cutting off the Fourier series at finite 
N, Eqs. (||) and (||) will not be satisfied identically, but 
instead there will appear two functions E\ and E2 such 
that 

A 5 V = Ei. # 0, B s i/> = E 2 ^0 

where E\ , E2 should be small in a sense that has yet to be 
determined. From the Fourier series expansion it can be 
seen, that increasing N will increase \E\\ too, since the 
conformal Laplacian annihilates every single term in the 
above Fourier series, L g <p n cos (wyz) — 0, for <p n known 
exactly. E\ is thus simply the sum of all numerical in- 
accuracies of the individual Fourier terms, and the sole 
purpose of increasing TV is to yield better satisfaction of 
the boundary conditions, that is to decrease E%. 

Let us deal with the smallness of E\ first. In the case 
of a vacuum spacetime, where the scalar curvature has to 
vanish, there does not exist an a priori measure for E\. 
Following jn| the error can be written as R — 16np res , 
where p res can be interpreted as the mass-energy density 
of matter for a non-vacuum spacetime. The criterion of 
accuracy derived from this viewpoint is the smallness of 
the total amount of mass represented by p res as compared 
to the ADM-mass. 

The residual of the Hamiltonian constraint is defined 
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FIG. 2. The ADM mass is plotted as a function of 7, the 
insert shows the pole of m at 7 = 70. The value for 7 = is 
marked by an x . 



D. Tests of Numerical Accuracy 

For every approximate solution 

N 

4>n = ^ ^n(r) cos {n-fz) 

n=0 



Pres = T7, lp L S ljj 

lbn 



yielding the residual mass 



I Pre 



\dV 



,\ip e d 3 x. 



Since the solution has been constructed in toroidal co- 
ordinates, one gets better results, if derivatives are com- 
puted with respect to these, so m res is rather evaluated 
in the form 



= — / ijjn- l ' 2 \L g G\fdfdzdip. 

2tT Jy 



The volume of integration has to be chosen carefully. 
Since the domain of computation is infinite with respect 
to the physical metric, integration over the whole com- 
putational domain would yield an infinite result, the in- 
tegration was thus restricted to \/r 2 + z 2 < 50 m. 

An alternative way to control the numerical error that 
is especially useful to control fulfillment of the boundary 
conditions can be constructed by putting the ADM mass 
into an explicitly positive form analogous to the expres- 
sion derived for Brill-waves in |). From Eq. (§) 

we get 

(d is the flat derivative operator) 



Q 



4> 



d 2 logVj + (<91og^) 2 



G 



which can be integrated with respect to the flat back- 
ground metric to give 



n a d a tp 



The boundary dM consists of spatial infinity and the 
boundary torus S, inserting the definition of the ADM 
mass to yields 



m 




n a d a j) 
1> 



dA 



M 



dtp 
I 7 



dV 



Using the boundary condition n a d a ip + ipp/4 — allows 
to calculate the first integral explicitly as 



1 

Z7T 



i> 8tt ' 



resulting in 



7T f r d^y __ 

777 = - + — / — dV =: m E + m v . 
2 27T J M V ^ / 



(31) 



The minimal mass of the metrics in the sequence consid- 
ered is thus bounded from below by m > n/2, the nu- 
merically calculated minimal value of the mass is 3.623 
at 7 = 0.628. 

In table (I) the relative errors of the numerically cal- 
culated expressions and m res /m and the choices of 
N are listed along with other results for some values of 
7- 

As the coordinate range of z approaches the entire real 
line for 7 — + 0, the cosine series (^7j) has to be replaced by 
a cosine integral that can be written down explicitly as 
is discussed in Appendix (jA]), accordingly the cutoff pa- 
rameter N of the Fourier series has to be increased for 7 
decreasing to keep the accuracy approximately constant. 
For practical reasons only few numerical computations 
were done for 7 < 0.15, so that the region < 7 < 0.15 
is left out of the figures displayed. 



where E is now taken to be our minimal toroidal bound- 
ary surface specified by f = 1 in the physical geom- 
etry, As is the Laplacian on E, R a b the Ricci cur- 
vature of g on E, and p a b the extrinsic curvature of 
E. Stability holds iff all eigenvalues of the operator 
— (Ax; + (R a bn a n b + p a bp ab )) are nonnegative for a min- 
imal surface. If E is of toroidal topology it is known Q| 
that it can only be stable if the minimal torus is totally 
geodesic, which it cannot be in our case as was concluded 
from the form of p a b given in Eq. ( |l4| ) . Nevertheless it 
is interesting to compute the spectrum (or at least the 
first few eigenvalues, since this has to be done numeri- 
cally). Using the Gauss law and the constraint conditions 
R = p = 0, and denoting the scalar curvature of E as 2 R 
the eigenvalue problem becomes 

- A s 0+i( 2 i?-;5 afc p Qb )<^7/0. 

Straightforward computation yields 



1 d 2 d 2 

1 — 7 2 dz 2 dip 2 



A s 


-G- 4 ( 




4G~ 4 


2 R 




1 - 7 2 


b p ab 





-2\-l 



With some rearrangement of terms the eigenvalue equa- 
tion becomes 



d 2 (j) 
' dz 2 



(1-7 2 ) 



9^ 

dip 2 



V(/> = i/G 4 (l - 7 2 



V=-- + 2 {G- 2 G,- Z 2 - G^Gw ) , 

which has to be solved with periodic boundary conditions 
in z and ip. Since the geometry is axially symmetric, so 
is G, and the equation separates. The eigenfunctions can 
be written as 



IV. INSTABILITY OF THE MINIMAL TORI 

An interesting question concerning surfaces of minimal 
area is whether they are stable - that is true minimal 
surfaces. Stability of a minimal surface E means that for 
all possible one-parameter (normal) variations E t of the 
surface E = Eo the second derivative of the area satisfies 
A"(0) > 0. A normal variation can be written as <f>n a 
where n a is a smooth unit normal vector field along Eo 
and <f> is a smooth function on E. The second variation 
of the area is then given by 



A"(0) 



[A E + (/UnV + p ab p ab 



)<j>]dA, 



4>(z, ip) = u nu (z)v n ((p), 

where v n (ip) = sin 77(1^9 + a) are the axial eigenfunctions 
degenerate in the axial phase shift a and for every n = 
0, 1, 2, . . . one gets an equation for the u nv (z): 



+ ((1 - 7 > 2 + V) u nu = Vn G\l - 7 2 )7 



(32) 



which is a regular self-adjoint Sturm-Liouville problem 
with periodic boundary conditions for 7 ^ 0,7q. The 
spectrum is thus purely discrete and bounded below (see 
e.g. Ref. 0, Sec. 1.4). 

As a resu t of using periodic boundary conditions the 
eigenfunctions are grouped in pairs by their number of 



7 



zeros (see e.g. Theorem 1.4.2 of Ref. jL7|): The eigen- 
functions Uk, k=0,l,2. . . , with eigenvalue i>k have k zeros 
for k even, and k + 1 zeros for k odd, that is u\ has no 
zero, U2 and 113 both have two zeros, U4 and 1*5 four, etc., 
where the cndpoints of the interval are identified. 

For practical purposes it is convenient to convert the 
periodic boundary conditions to unmixed boundary con- 
ditions, which is made possible by noting that for Eq.(^) 
it is sufficient to consider only eigenfunctions that are 
with respect to z = 0. If a periodic solution u exists 
that is neither symmetric nor antisymmetric, then the 
whole solution space is spanned by the linear combina- 
tions u(x) ± u{—x) and the general solution is periodic, 
since with u(x) being a periodic solution, so is u{— x). 
The eigenvalue is of multiplicity two in these cases, which 
show up in figures (^,^) as the crossing points of branches 
associated with symmetric and antisymmetric eigenfunc- 
tions that have both the same number of zeros. 

The problem of finding solutions with periodic bound- 
ary conditions can thus be split into the easier problems 
of finding all periodic symmetric and antisymmetric so- 
lutions, equivalent to imposing the unmixed boundary 
conditions 



u'(0) = u' Q 



= for symmetric solutions, 



u(0) = u I — I = for antisymmetric solutions. 
V7/ 

so that standard methods for Sturm-Liouville problems 
can be applied. Here routine d02kef from the NAG For- 
tran Library was chosen to calculate eigenvalues and 
eigenfunctions, as was done already for Eq. (p5|). 

The results are given in figures (||,[|) for n = 0, 1. For 
all values of 7 one can find at least one negative eigen- 
value, demonstrating that the tori are not stable, and are 
thus only saddle points of the area functional. In particu- 
lar the lowest eigenvalue for n — is found to be negative 
for all 7, and all other eigenvalues for n = become neg- 
ative as well for 7 small enough, see Fig. (^). For n = 1 
only the lowest eigenvalue is negative for 7 > 0.71, see 
Fig. For n = 2 the situation is similar: again only 
the lowest eigenvalue is negative, but in the much smaller 
region 7 > 0.97. For n > 2 the situation is not totally 
clear: if negative eigenvalues appear, they have to occur 
very close to 70. By explicit computation of the area 
of neighboring tori with f = const., it was found that 
the minimal tori have minimal area for 7 < 0.891, and 
maximal area for 7 > 0.891 among neighboring tori of 
constant f. 

In the limiting case 7 — 70 the Green function G di- 
verges, and all eigenvalues v go to zero. In the other 
limiting case, 7 = 0, the coordinate range is the entire 
real line, and the Sturm-Liouville problem becomes sin- 
gular. The potential V falls off exponentially to the value 
n 2 — 1/4 for large \z\ - note that this value is negative 
for n — and positive for n > 0. While nothing spe- 
cial happens to the part of the spectrum with n > 0, see 



Fig. (0) for n — 1, the situation is different for n = 0. 
The general solution u to Eq. ( |32| ) for 7 = n = and 
arbitrary v is asymptotically trigonometric (see e.g. Ref. 
ifLsf , section 5.11), that is 

u = asin(z + f3) + 0(l/z), 

which is oscillatory with finite amplitude near z = ±00. 
Periodic boundary conditions thus can not be imposed on 
these solutions, for small 7 all eigenvalues decrease un- 
boundedly, marking a growing instability of the minimal 
surface with respect to axially symmetric deformations. 
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FIG. 3. The first five eigenvalues vqi (n — 0) are plot- 
ted as functions of 7 from top down (I = 1,...,5). The 
curves labelled 1,2,3 are the first three branches of eigen- 
values with symmetric eigenfunctions, while l',2' label the 
first two branches of eigenvalues with antisymmetric eigen- 
functions. 



V. APPARENT HORIZONS 

As was shown before, the minimal tori considered here 
are all unstable, so that they locally can be deformed to 
an outer trapped surface [Q . By asymptotic flatness the 
region of compact outer trapped surfaces has to be finite, 
the outer boundary is called apparent horizon, which is 
a stable marginally outer trapped surface (see Ref. 0] 
and references cited therein). Noting that g is analytic 
by construction of the conformal factor as the solution of 
Eqs. ( p8f]30| ) and not flat (by analyticity it would have to 
be flat everywhere if it were flat in some region, but this 
can not be, since the mass is positive by Eq. (|3l])), it can 
be concluded from Ref. |Q that the apparent horizon has 
spherical topology. 

The existence and properties of apparent horizons have 
received some interest in the context of the cosmic cen- 



sorship hypothesis |19|, where they play a double role as 
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FIG. 4. The first five eigenvalues V\i in — 1) are plotted 
as functions of 7 from top down (I — 1, ... ,5), the values for 
7 = are marked by x's. The curves labelled 1,2,3 are the 
first three branches of eigenvalues with symmetric eigenfunc- 
tions, while l',2' label the first two branches of eigenvalues 
with antisymmetric eigenfunctions. 



in Penrose's singularity theorem (To| and also as a local 
(albeit only weak) analogue of the global concept of an 
event horizon. The second is practically exploited in the 
apparent horizon boundary condition || used in numer- 
ical relativity and also was used in a new approach for 
the definition of a black hole itself [^(J . 

The apparent horizons were located numerically and 
analyzed with respect to the Penrose and hoop conjec- 
tures. Penrose |li]] formulated the conjecture that if a 
collapsing body satisfies the inequality 



m < 



A 

16^' 



where m is the ADM mass and A the area of a trapped 
surface, then cosmic censorship should be violated. As 
pointed out by Horowitz if A is the minimal area 
that surrounds a trapped surface (in the time symmetric 
case identical to the apparent horizon) and satisfies the 
above inequality, then cosmic censorship does not hold. 
For the apparent horizons found here the inequality is 
always violated, the Penrose conjecture thus confirmed, 
see figure (||). 

The hoop conjecture, due to Thorne Jll]], states that 
an apparent horizon will exist if and only if a mass m gets 
compacted into a region such that the circumference C 
satisfies C < Airm in every direction. The exact defini- 
tion of what should be meant by m and C is left open. 
For a precise definition of m one would need a notion of 
quasilocal mass which is available only for the spherically 
symmetric case (no gravitational waves), here the ADM 
mass was chosen. Also for a general definition of the 



circumference various definitions would be possible. In 
an axially symmetric geometry which is also symmetric 
with respect to reflection at the equator, straightforward 
definitions of circumference exist along the equator (C e ) 
and along a meridian circle (C p ): 



ri 2 , 

C p = 4 dd^VR' 2 + R 2 , 
Jo 

where R is the radial coordinate in standard polar coor- 
dinates (R,9,ip). The results for are given in Fig. 
(|) and table (I). 

For the problem of locating apparent horizons in the 
axially symmetric case a variety of different methods is 
available by now. If one considers configurations which 
are symmetric (or almost symmetric) with respect to re- 
flection at the equator, one can reasonably hope that the 
apparent horizon will be expressible as a function R(8). 
This gives the ODE 



b" Rl ' 2 



4R 2 tH+Ri f4Ucot, 



2R = 0, 



which has to be solved with the boundary conditions 

R'(0) = R'(tt/2) = 0. 

This is done numerically using a shooting and matching 
method (NAG Jl5| routines d02agf and d02ebf were used 
in two different codes for horizon finding). 

The results are shown in figures (|), (§) and table (I). 
For 7 — * 70 the apparent horizon becomes spherical and 
pinches off to infinity as has been discussed for the anal- 
ogous case when the Yamabe number goes to zero for a 
manifold without boundary in Ref. ||. For 7 < 70 the 
apparent horizon is oblate, but not very strongly. 

VI. DISCUSSION 

In a previous paper [ fl2"[ the notion of toroidal confor- 
mal symmetry was introduced to label a class of initial 
data where a U(l) x U(l) conformal isometry is present. 
Working on a compactified background metric allows to 
easily exploit this conformal symmetry to decompose the 
Lichnerowicz conformal factor in a double Fourier series 
on the group orbits and give the solution in terms of a 
countable family of uncoupled ODEs on the orbit space. 
The present paper applies these methods to construct 
and analyze initial data containing a marginally outer 
trapped torus. Special emphasis is given to the investiga- 
tion of the limiting cases when the sequence of solutions 
intersects the boundary of conformal superspace. 

To get a better understanding of the results it seems 
useful to look at the usually studied case of vacuum 
asymptotically flat data without boundary, where the 
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FIG. 5. Check of the Penrose isoperimetric inequality, the FIG. 6. Equatorial and polar hoop quotients, the values for 
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cases understood best are the time-symmetric, notably 
the axially symmetric subclass known as Brill- waves j22| . 

In order to simplify this comparison I will first discuss 
generalizations of the data considered in this paper. The 
obvious generalization is to drop conformal flatness and 
consider deformations of the base metric g of the form 



= e 2A 9 (r) 



dr 



+ {l-i z r 2 )d? I +f 2 dip 2 , (33) 



which is the analogue of the Brill-wave ansatz [Q while 
keeping the conformal symmetry, and is currently inves- 
tigated. The general axially symmetric case corresponds 
to form functions q(f, z) that depend on both f and z. 
Within the present ansatz one could consider certain non- 
axially symmetric situations by placing A off the axes 
(the sets of fixed points of the isometries), the physical 
metric then has no isometries whatsoever - since an isom- 
etry of g carries over to g only if A is a fixed point of this 
isometry p3|. 

The usually studied case is recovered, finally, by impos- 
ing regularity on all of R 3 instead of the p — boundary 
condition, this case has been discussed in jl^], the pub- 
lication of numerical results is in preparation. 

For these traditionally considered data sequences of 
metrics which approach the boundary of the space of ini- 
tial data have been studied both analytically and numer- 
ically and tested for apparent horizons. Two comple- 
mentary paths have been taken, an overview is given in 
p2| . Sequences of metrics for which the Yamabe num- 
ber goes to zero have been studied analytically by Beig 
and 6 Murchadha [§, who proved that asymptotically 
the mass diverges in indirect proportion to the Yamabe 
number and that an apparent horizon would always ap- 
pear before the mass diverges. Abrahams et al. p3] on 



the other hand studied families of metrics of fixed mass 
with a singular limit for which no apparent horizon ap- 
peared even in the limiting case. 

The limiting case 7 — > 70 of the present paper is sim- 
ilar to the one where the Yamabe number goes to zero. 
As the Sobolev quotient defined in Eq. (|ll|) is just an 
extension of the Yamabe number including an additional 
boundary term, the mathematical and physical situation 
is analogous. As the mass diverges, the apparent hori- 
zon becomes more and more spherical and pinches off to 
infinity. Going beyond that critical point the mass be- 
comes infinitely negative. Whereas in the no-boundary 
case the amplitude parameter A of Eq. (|33| ) can be in- 
creased without bound, which gives rise to an infinite 
number of analogous singularities in the mass function 
(to be discussed in forthcoming work), here the param- 
eter 7 is bounded from above by 7 < 1, and the mass 
appears to stay negative for 7 > 70- A proof that the 
Sobolev quotient will become negative was unsuccessful 
due to the finite range of 7: In the Brill- wave case an 
estimate can be obtained with simple test-functions by 
just increasing the amplitude enough [^2| - here the same 
trick can not be carried out. 70 comes in fact very close 
to the upper boundary 7=1, and so a very good choice 
of test function would have to be made. 

The other case where the sequence of solutions hits 
the boundary of conformal superspace, 7 — > 0, is fun- 
damentally different. While in the first case the base 
geometry stays regular, the situation when the boundary 
torus touches the axis can be considered as singular, nev- 
ertheless the mass and location of the apparent horizon 
remain finite, but a curvature singularity develops at the 
origin. 

Table (I) lists some numerical results for different val- 
ues of 7. 
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A0 = 0, 



TABLE I. Results from numerical calculations, numbers 
are rounded to the last digit. 



The technique of using conformal symmetry together 
with compactification had several technical advantages. 
First of all, the ODEs ( p8| ) and (^5|) resulting from the 
Hamiltonian constraint and existence question can be 
handled with relatively small numerical expertise and 
desktop workstations. The use of the NAG numerical 
Fortran library p"5[ yields an excellent documentation of 
all numerical routines that have been used. 

Another point in favor of the method is that it works 
well for strong data - that is near 7 = 70 in our sequence, 
only the n = term in the Fourier series blows up when 
the mass does. 

If there are additional asymptotic regions (representing 
black holes), discrete symmetries result from an arrange- 
ment compatible to the conformal isometry as discussed 

in @- 

It is hoped, that the technical approach taken here 



and in [ 
mal Kil 



12[ , namely to use compactification and confor- 
ing vector fields may also prove useful in other 



situations. 
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APPENDIX A: THE EXPLICIT LIMITING 
SOLUTION 

In the limiting situation 7 = the constraint equation 
(B2n can be solved explicitly. Making the ansatz 



G 



v/2 



\/f 2 + z 2 



(analogous to (p6h), we get 



= (l + 2 2 )- 3 / 2 -^(l + Z 2 



-1/2 



Or 



from Eqs. (f22|) and 

We solve by cosine transformation: 



(r,z) = - 

7T 



cos(uJz)f u] {r) dui, 



resulting in a family of Bessel equations for regular func- 
tions f u {f) parameterized by u>: 



The regular solutions are 

U{r) = c(uj)I (ujf), 

where /„ and K n used below are modified Bessel func- 
tions of order n and c(uj) is determined from the cosine 
transformed boundary condition 



cos(luz) 

((i + z 2 r 3/2 -^(i + - 2 r 1/2 ) ** 



2 4\coK 1 (to) - K (lo) 



_ 1 

~ 4 

as 

c{lo) , 

7T AloI\{lo) + Io{ui) 
We get an explicit formula for the Green function, 

OTA * ■ * 



\/r 2 + z 2 



\Jf 2 + z 1 



• . 4u{Ki{u)) -K (u) j , — \ 7 , A1 , 
«»M , r 1 \ ~Z t ( \ L h{ur)dw. Al 



<f> is a regular function of r, z, the integrand diverges 
only logarithmically for small u> and exhibits exponential 
falloff for large \z\. Its value at the point at infinity gives 
the ADM mass m = 3.877 by m = 2V2</>(A). 

Alternatively the Green function can be written as 



p(u>,r) = Auj(Ii(uj)K (ujr) + Ki(uj)I (ujr)) 
+K (uf)I (u) - I (wf)XoM, 
q(ui) = 4uli(u) + Iq(u). 



using the identity 
1 



y/r 2 + 



cos(w z)Ko(ujr) du>. 
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Positivity of G(f, z) was checked by numerical evalua- 
tion, the behavior of G for large \z\ can be seen analyti- 
cally by integration of (A2) in the complex plane. Both 
p and q are positive and symmetric on the real line, but 
q has infinitely many zeros for purely imaginary argu- 
ments, leading to corresponding poles of p/q. Observ- 
ing, that the factor exp (— Im(w)z) makes contributions 
to the integral small for large z in the upper complex half 
plane (for z — ► — oo the lower half respectively) yields an 
asymptotic formula for the integral by choosing a path 
that encircles only the first pole (an asymptotic series 
results from summing over all poles), 

G(f, z) — c(f) exp (— a\z\), 

where a = 0.6856 is the first zero of q(to) on the imagi- 
nary axis and 



c(f) = 4 i Residual ( ,u = 0.6856 



Some information is also available over WWW, at 



http : //www .nag . co.uk/ 
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